Weak interaction rates for Kr and Sr waiting-point nuclei under rp-process conditions 
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Weak interaction rates are studied in neutron deficient Kr and Sr waiting-point isotopes in ranges 
of densities and temperatures relevant for the rp process. The nuclear structure is described within 
a microscopic model (deformed QRPA) that reproduces not only the half-lives but also the Gamow- 
Teller strength distributions recently measured. The various sensitivities of the decay rates to 
both density and temperature are discussed. Continuum electron capture is shown to contribute 
significantly to the weak rates at rp-process conditions. 
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INTRODUCTION 



' Nuclear physics is an essential piece in the present understanding of many astrophysical processes related to the 
energy generation and nucleosynthesis in stellar scenarios. The properties of both stable and exotic nuclei play different 
\^ • roles at different stages of stellar evolution. While the former are relevant to describe steady regimes, the latter are 
implicated in the description of explosive events. Network calculations and astrophysical models rely on the quality of 
the input information, in particular, on the quality of the nuclear input. Unfortunately, the experimental information 
available for exotic nuclei is still very short and thus, most of the astrophysical simulations of violent phenomena 
must be based on nuclear model predictions of limited quality and accuracy. This is the case of the X-ray bursts, 
O ' generated by a thermonuclear runaway in the hydrogen-rich environment of an accreting compact object (typically a 
^ , neutron star) which is fed from a binary companion (typically a red giant). The ignition takes place on the surface 
of the neutron star at high densities (p) and temperatures (T), and eventually peak conditions of T =1-3 GK and 
p — 10^-10'' g cm~^ are reached [l|, The mechanism leading to X-ray bursts is the rapid proton (rp) capture 
process [J [1, S] , which is characterized by proton capture reaction rates that are orders of magnitude faster than any 
other competing process, in particular /9-decay. It produces rapid nucleosynthesis on the proton-rich side of stability 
toward heavier proton-rich nuclei and the energy observed in X-ray bursts. 
■ Nuclear reaction network calculations (the set of differential equations for the various isotopic abundances) are 
performed [J 0) S 0] to follow the time evolution of the isotopic abundances, to determine the amount of energy 
released by nuclear reactions, and to find the reaction path for the rp process. In general, the reaction path follows a 
series of fast proton-capture reactions until further proton capture is inhibited. Then the reaction flow has to wait for 
a relatively slow /3-decay and the respective nucleus is called a waiting point (WP). The half- lives of the WP nuclei 
\ along the reaction path determine the time scale of the nucleosynthesis process and the produced isotopic abundances. 
• • I In this respect, the weak decay rates of neutron-deficient medium-mass nuclei under stellar conditions play a relevant 
role to understand the rp process. The nuclear structure problem involved in the calculation of these rates must 
be treated in a reliable way. In particular, this implies that the nuclear models should be able to describe at least 
5^ the experimental information available on the decay properties (Gamow- Teller strength distributions and /3-decay 
5^ , half-lives) measured under terrestrial conditions. Although these decay properties may be different at the high p 
and T existing in rp-process scenarios, success in their description under terrestrial conditions is a requirement for 
a reliable calculation of the weak decay rates in more general conditions. With this aim in mind, we study in this 
letter the dependence of the decay rates on both p and T , using a nuclear model that has been tested successfully 
and reproduces the experimental information available on both bulk and decay properties of medium-mass nuclei very 
reasonably. This model is the quasiparticle random phase approximation (QRPA). Here, we focus our attention to the 
deformed WP isotopes '^2,74-[^j. 76,78gj.^ where the Gamow- Teller (GT) strength distributions have been measured 
with high accuracy by /3-decay over most of the Q- window [5|, (a, III and are well reproduced by theoretical calculations 
based on deformed QRPA ^8j|. While the half-lives give only a limited information of the decay (different strength 
distributions may lead to the same half-life) , the strength distribution contains all the information. It is the first time 
that the weak decay rates under stellar rp conditions can be studied using a nuclear structure model that reproduces 
the GT strength distributions and half- lives under terrestrial conditions [1, Q . 
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There are several distinctions between terrestrial and stellar decay rates caused by the effect of high p and T. 
The main effect of T is related to the thermal population of excited states in the decaying nucleus, accompanied by 
the corresponding depopulation of the ground states. It has been shown [l^ that the weak decay rates of excited 
states can be significantly different from those of the ground state, enhancing the total decay rates. A case by case 
consideration is needed. Another effect related to the high p and T comes from the fact that atoms in these scenarios 
will be completely ionized and the electrons will be no longer bound to the nuclei, but forming a degenerate plasma 
obeying a Fermi-Dirac distribution. This opens the possibility for continuum electron captures (cEC). These effects 
make weak interaction rates in the stellar interior sensitive functions of T and p, with T = 1.5 GK and p = 10^ g 
cm~'^, as the most relevant conditions for the rp process 



II. WEAK DECAY RATES IN STELLAR SCENARIOS 



The general formalism to calculate weak interaction rates in stellar environments as a function of p and T was 
introduced in the pioneering work of Fuller, Fowler, and Newman pT| . Further improvements have come mainly on 
the nuclear structure aspect, either from the Shell Model [13; [3 or from QRPA [141] . 

The decay rate of the parent nucleus is given by [ll| 

A = ^ A, ^:^e-^-/('=^) ; G = ^(2J, + l)e~^^/('^^-> , (1) 

i i 

where Jj {Ei) is the angular momentum (excitation energy) of the parent nucleus state i, and thermal equilibrium is 
assumed. In principle, the sum extends over all populated states in the parent nucleus up to the proton separation 
energy. However, since the range of T for the rp-process peaks at T = 1.5 GK {kT ~ 300 keV), only a few low-lying 
excited states are expected to contribute significantly in the decay. Specifically, we consider in this work all the 
(collective) low-lying excited states below 1 MeV [l^[3| (exception made of the £^4+ = 782 keV in ^^Sr), Eq+ = 671 
keV and E2+ = 709 keV in ^^Kr; Eo+ = 509 keV and E2+ = 456 keV in ^^Kr; E2+ = 261 keV in ^^Sr; and E2+ = 279 
keV in ^^Sr. Then, all these states are considered to be thermally populated although for T = 1.5 GK the population 
is negligible beyond 500 keV. In Fig. [T] one can see the population coefficients in Eq. ([1]) as a function of the excitation 
energy (left) and temperature (right). As an example, the population of the lowest of these states (£'2+ — 261 keV in 
^^Sr) at r=1.5 GK is about 12%. 
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FIG. 1: Population coefficients as a function of the excitation energy (left) and temperature (right). The energies of the 
low-lying excited states of the isotopes considered in this work are shown explicitly. 

The decay rate for the parent state i is given by 
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where the sum extends over ah the states in the final nucleus reached in the decaying process. The rate Xif from the 
initial state i to the final state / is given by 

= (3) 

where D = 6146 s. This expression is decomposed into a nuclear structure part Bif and a phase space factor ^if. 
Bif contains the transition probabilities for allowed Fermi and GT transitions, 

B.,f = B.,f{GT)+B,j{F), (4) 

with 
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cff 



{9A/gv)cS = 0.74(5yi/(7y)barc is an effective quenched value. 

B.f{F) = ^^(/ll ^ ti\\^r = t{t + 1) - t,,,Uf , (6) 

k 

where t is the isospin, which is assumed to be a good quantum number. For /?''" transitions and very neutron 
deficient nuclei with Z > N, the superallowed Fermi transition to the isobaric analog state has to be taken into 
account. However, for N > Z nuclei, Fermi transitions are not possible except for isospin impurities. The discussion 
is therefore limited to the decay via allowed GT transitions. 

The theoretical formalism used here to calculate the GT strengths is based on th e Q RPA. The quasiparticle basis 
corresponds to a deformed selfconsistent Skyrme Hartree-Fock calculation with SLy4 [i7| force and pairing correlations 
treated in BCS approximation. The residual interactions are spin-isospin forces in both particle-hole and particle- 
particle channels. The initial and final states in the laboratory frame are expressed in terms of the intrinsic states 
using the Bohr-Mottelson factorization which is a very good approximation for well deformed nuclei. The effect 
of all these ingredients (deformation, pairing, residual forces) on the decay rates has been studied elsewhere 11911 
and their relevance in this mass region has been stressed. More details of the formalism can be found in Rcf. [2(j |. 
The present formalism has been shown to provide a good description of the decay properties of nuclei in the mass 
region ^ ^ 70. In particular, the parameters for residual forces used here are the same used in Ref. Q, where good 
agreement was obtained with the experimental GT strength distributions and /3-decay half-lives. The success of this 
theoretical formalism in reproducing terrestrial decay rates supports its application to the calculation of the weak 
interaction rates in stellar matter. 

The P'^ and cEC phase space integrals are given by 



where 



<i>.; = $^;^ + <i>f;, (7) 

= r uj^u:^-l{Q,f+ujfF{Z,u:) 

J tjjf 

x5eM[l-5,(g,/-fa;)]dw, (8) 

for continuum electron capture, and 

= / oj^uj-^ - l{Q,f - ujfF{-Z + l,cj) 

^[l~Sp{uj)][l-S,{Q^f-Lo)]dw, (9) 

for positron emission. In these expressions uj is the total energy of the electron (positron) in units of TOgC^ and F{Z, uj) 
is the Fermi function [2l| that takes into account the distortion of the /3-particle wave function due to the Coulomb 
interaction. 
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FIG. 2: Electron distributions as a function of the total electron energy for various T at fixed p (left panel) and for various p 
at fixed T (right panel). 



Qif = 



Ta„c' 



(10) 



is the total energy available in the decay in units of m^c? . It is written in terms of the nuclear masses of parent 
(Mp) and daughter {Ma) nuclei and their excitation energies Ei and Ef^ respectively. In the cEC factor, the lower 
integration limit is given by = 1 if Qi/ > —1, or = \Qif \ H Qif < —1. ^e, Sp, and S,^, are the electron, positron, 
and neutrino distribution functions, respectively. Its presence inhibits or enhances the phase space available. In rp 
scenarios the commonly accepted assumptions [ll, [ll|, IH, [H, state that Si, = 0, since neutrinos can escape freely 
from the interior of the star without blocking their emission in the capture or decay processes. Positron distributions 
become only important at higher T {kT > 1 MeV) when they appear via pair creation. At the temperatures considered 
here we take Sp — 0. The electron distribution is described by a Fermi-Dirac distribution 
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FIG. 3: Phase space factors for Kr as a function of the total energy available in the decay for various T (GK) and p (g cm~^). 
Electron capture factors are calculated from orbital electrons in the case T — (a.) and from continuum electrons in the case 
of nonzero T (b,c). 



Se 



1 



exp [{lo ~ /{kT)] + 1 



(11) 
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assuming that nuclei at these temperatures are fully ionized and the electrons are not bound to nuclei. The chemical 
potentials /ig are determined from p and T. Tables for them can be found in Ref. [ll[ for selected choices of T and p. 

Fig. [2] shows Se{bj) for various temperatures and electron densities pY^ [p is the baryon density and Ye stands for 
the ratio of the electron to the baryon numbers) around the peak rp conditions. One can see in Fig.[^that Se becomes 
more important at higher densities, as well as the relevance of electrons with high energies as T increases. These 
features establish the relative importance of the cEC rates as p and T change. 

Fig. [3] contains the phase space factors for /?+ decay and electron capture as a function of the total energy available 
Qif. In (a) one can see the usual phase factors at zero T, that is the_5+ in Eq. ([§]), which is independent of T and 



p and the orbital electron capture calculated as described in Ref. [2l|. In (b) and (c) one can also see the cEC 
factors as a function of p and T. In general, they increase with Qif and thus the decay rates are more sensitive to 
the strength Bif located at low excitation energies of the daughter nucleus {Ef). It is also interesting to notice the 
relative importance of both /J"*" decay and electron capture. The latter is always dominant at sufficiently low Qif, or 
correspondingly at sufhciently high excitation energies Ef. In this work the effects of bound-state electron capture 
are not included since they are insignificant over most of the ranges of T and p considered. 




T (GK) T (GK) T (GK) 



FIG. 4: Weak decay rates A (s ^) for ^^Kr as a function of T for various p. 
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FIG. 5: Same as in Fig. H but for '^^Sr. 



III. RESULTS 



Fig. [5] shows the decay rates versus T on the example of ""^Kr. In panel (a) one can see the decomposition of 
the total rate into their contributions from the decay of the the ground state 0+ — > 1+ and from the decay of the 



6 



excited states 0^ 1+ and 2+ l+,2+,3"'", which are neghgible except at very high T. Panel (b) contains the 
decomposition of the rates into their /3+ and cEC components evaluated at different densities. One can see that for 
densities larger than log]^Q(p) ^ 6.5, cEC becomes dominant at any T. For lower densities, rates are larger than 
cEC at low T, but lower at sufficiently high T. Since the /3+ rate is practically independent of both p and T, the total 
rate in panel (c) is determined by /3+ at low T and p, and by cEC otherwise. The gray area is the relevant range 
r=l-3 GK for the rp process. Fig. [5] shows the same rates as in Fig. HI but for ^^Sr. The only difference worth to 
mention with respect to ^^Kr is the relatively more important contributions of the 2+ excited state, as seen in panel 
(a), which is the result of a larger population of this state due to its lower excitation energy. 

Finally, Fig. [6] contains the half- lives (T1/2 = In 2/ A), including /3+ and cEC contributions, as a function of T at 
a fixed p = 10^ for the various isotopes. One can see the decrease of T1/2 as T increases. The decrease starts to 
be significant beyond T = 3 GK, which is outside of the relevant temperature range for the rp process. One should 
notice that the orbital electron capture has been ignored in this work in favor of continuum electron capture, but 
both should be considered to obtain a smooth transition toward terrestrial conditions. The calculation at T = with 
orbital electron capture contributions can be found in Ref. [sj , where comparison with experiment is also shown. 

1 — I — I — I I I U 
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FIG. 6: Half-lives as a function of T at a density p = lO'' g cm~^ for Kr and Sr isotopes. 



IV. CONCLUSIONS 

In conclusion, the weak decay rates of Kr and Sr WP nuclei have been investigated at temperatures and densities 
where the rp process takes place. It has been the first time that stellar weak rates have been calculated in this 
mass region with a nuclear structure model that reproduces both the half-lives and the more demanding GT strength 
distributions measured with high accuracy from /3-decay experiments. In this letter we have analyzed the relevant 
ingredients to describe the rates in a reliable way. We have studied: 1) The contributions to the decay rate coming 
from excited states in the parent nucleus which are populated as T raises. All the relevant states below 1 MeV 
have been included. It is found that below T = 3 GK their effect can be neglected, and thus, the decay from the 
ground state is already a good approximation for rp processes. This conclusion is in agreement with previous studies 
EE El. Nevertheless, one should pay special attention to the cases where the 2+ excited states are particularly 
low in energy because their contributions can be competitive at these temperatures; 2) The effect of the continuum 
electron capture rates. It is found that cEC rates are enhanced as T and p increase and they become comparable to 
the (3'^ decay rates at rp peak conditions. At slightly larger values of T and p, cEC dominates over decay. This 
point is important since cEC contributions have been neglected [l| in earlier evaluations of weak decay rates at rp 
conditions. 
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